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Abstract 

The thermal physics of a massless scalar field with a (j)'^ interaction is studied within 
screened perturbation theory (SPT). In this method the perturbative expansion is 
reorganized by adding and subtracting a mass term in the lagrangian. We con- 
sider several different mass prescriptions that generalize the one-loop gap equation 
to two-loop order. We calculate the pressure and entropy to three-loop order and the 
screening mass to two-loop order. In contrast to the weak-coupling expansion, the 
SPT-improved approximations appear to converge even for rather large values of the 
coupling constant. 
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I. INTRODUCTION 



If we have a weakly-coupled quantum field theory in equilibrium at temperature T, we 
should be able to use perturbation theory as a quantitative tool to study its properties. In 
the case of a massless theory with coupling constant g, the naive perturbative expansion 
in powers of breaks down because of collective effects such as screening. However, the 
perturbative expansion can be reorganized into a weak-coupling expansion in powers of g 
either by using resummation methods or alternatively by using effective field theory. It 
is reasonable to assume that this weak-coupling expansion provides a useful asymptotic 
expansion for sufficiently small values of g. 

Only in recent years has the calculational technology of thermal quantum field theory 
advanced to the point where this assumption can be tested. Unfortunately, the assumption 
seems to be false. One would expect the thermodynamic functions, such as the pressure, 
to be among the quantities with the best-behaved weak-coupling expansion, since collective 
effects are suppressed by several powers of g. However, in recent years, the thermodynamic 
functions have been calculated to order g^ for massless scalar theories abelian gauge 

theories 0,111, and nonabelian gauge theories The weak-coupling expansions show 

no sign of converging even for extremely small values of g. There is already a hint of the 
problem in the g^ correction, which has the opposite sign and is relatively large compared 
to the g^ coefficient. The large size of the g^ contribution is not necessarily fatal, since 
it is the first term that takes into account collective effects. An optimist might still hope 
that higher-order corrections would be well-behaved. This optimism has been dashed by the 
explicit calculation of the g'^ and g^ terms. 

For a massless scalar field theory with a (7^0*^/4! interaction, the weak-coupling expansion 
for the pressure to order g^ is pl- 



V = V; 



Meal l~\a + ^«^/^ + ^ (log ^ + 0.40) a' 
^^"^ ^log ^ - ^ " - 0-^2) a'/' + 0{a' log a) 



(1) 



where Pideai = (vr^/90)T'' is the pressure of an ideal gas of free massless bosons, a = 
g"^ (fi) / IGn"^ , and g{fi) is the MS coupling constant at the renormalization scale /i. In Fig. |1|, 
we show the successive perturbative approximations to P/Pideai as a function of g{2'7iT). 
Each partial sum is shown as a band obtained by varying fi from vrT to 47rT. To express 
g{n) in terms of g{2'KT), we use the numerical solution to the renormalization group equation 
/i^a = (3{a) with a five-loop beta function IQ]: 

d 1 7 

a a = 3a^ + 32.54a^ - 271.6a^ + 2848.6a^ . (2) 

Ofi 3 

The lack of convergence of the weak-coupling expansion is evident in Fig. |1|. The band 
obtained by varying /i by a factor of two is not necessarily a good measure of the error, but 
it is certainly a lower bound on the theoretical error. Another indicator of the theoretical 
error is the deviation between successive approximations. We can infer from Fig. [l| that the 
error grows rapidly when g{27iT) exceeds 1.5. 
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FIG. 1. Weak-coupling expansion to orders 5^, , g^, and g^ for the pressure normalized to 
that of an ideal gas as a function of g{2TrT). 



A similar behavior can be seen in the weak-coupling expansion for the screening mass, 
which has been calculated to next-to- next-to-leading order in g [H]: 



aT^l 1 - V6a'/' - 



3 log 



2nT 



2 log a -6.4341 



a + 0{a 



3/2 N 



(3) 



In Fig. 0, we show the screening mass normalized to the leading order result m^o = 
g{2T{T)T / \/2A as a function of g(27rT), for each of the three successive approximations to 
m^. The bands correspond to varying fi from irT to AttT. The poor convergence is again 
evident. The pattern is similar to that in Fig. |I|, with a large deviation between the order-^f^ 
and order-(7^ approximations and a large increase in the size of the band for g^. 

There are many possibilities for reorganizing the weak-coupling expansion to improve 
its convergence. One possibility is to use Fade approximants 0. This method is limited 
to observables like the pressure, for which several terms in the weak-coupling expansion 
are known. Its application is further complicated by the appearance of logarithms of the 
coupling constant in the coefficients of the weak-coupling expansion. However, the greatest 
problem with Fade approximants is that, with no understanding of the analytic behavior of 
V at strong coupling, it is little more than a numerological recipe. 

m 



Fer- 



An alternative with greater physical motivation is a self-consistent approach 
turbation theory can be reorganized by expressing the free energy as a stationary point of a 
functional fl of the exact self-energy function n(j9o, p) called the thermodynamic potential 



T^. Since the exact self-energy is not known, 11 can be regarded as a variational function. 
The "^-derivable" prescription of Baym [jlO[ is to truncate the perturbative expansion for 
the thermodynamic potential Q and to determine 11 self-consistently as a stationary point 
of Q. This gives an integral equation for 11 which is difficult to solve numerically, except in 
cases where 11 is momentum independent. In relativistic field theories, there are additional 
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FIG. 2. Weak-coupling expansion to orders , g^, and for the screening mass normalized to 
the leading-order expression as a function of g{2TrT). 



complications from ultraviolet divergences. A more tractable approach is to find an approx- 
imate solution to the integral equations that is accurate only in the weak-coupling limit. 
Such an approach has been applied by Blaizot, lancu, and Rebhan to massless scalar field 
theories and gauge theories fT^JTSH . 

Another approach that is also variational in spirit is screened perturbation theory (SPT) 
introduced by Karsch, Patkos and Petreczky [T^. This approach is less ambitious than the 
^-derivable approach. Instead of introducing a variational function, it introduces a single 
variational parameter m. This parameter has a simple and obvious physical interpretation 
as a thermal mass. The advantage of screened perturbation theory is that it is very easy 
to apply. Higher order corrections are tractable, so one can test whether it improves the 
convergence of the weak-coupling expansion. Karsch, Patkos and Petreczky applied screened 
perturbation theory to a massless scalar field theory with a 0^ interaction, computing the 
two- loop pressure and the three-loop pressure in the large- limit. In both cases, they used 
a one-loop gap equation as their prescription for the mass. Their three-loop calculation was 
not a very stringent test of the method, because the large-A^ limit suppresses self-energy 
diagrams that depend on the momentum. 

In this paper, we present a thorough study of screened perturbation theory for a massless 
scalar field theory with a (p^ interaction. We calculate the pressure and entropy to three 
loops and the screening mass to two loops using SPT. We consider several generalizations 
of the one-loop gap equation to two loops. Inserting the solutions to the gap equations for 
m into the SPT expansions, we obtain the SPT-improved approximations to the pressure, 
the screening mass, and the entropy. 

The paper is organized as follows. In section ||, we describe the systematics of screened 
perturbation theory. In section we discuss the possible prescriptions that can be used to 
generalize the one-loop gap equation to higher orders. We calculate the free energy to three- 
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loop order in section |I^ and the screening mass to two- loop order in section 0. In section 0, 
we study three generalizations of the one- loop gap equation to two- loop order. In section VII , 
we study the convergence of the SPT-improved results for the pressure, screening mass, and 
entropy. In section [VIII| , we summarize and conclude. We have collected the necessary 
sum-integrals in an appendix. 



II. SCREENED PERTURBATION THEORY 

The lagrangian density for a massless scalar field with a (f)'^ interaction is 

£ = ^d,<j>dy - ^gY + A£ , (4) 

where g is the coupling constant and AC includes counterterms. The conventional pertur- 
bative expansion in powers of g'^ generates ultraviolet divergences, and the counterterm AL 
must be adjusted to cancel the divergences order by order in g"^. If we use dimensional regu- 
larization in (i = 3 — 2e spatial dimensions and minimal subtraction to remove the ultraviolet 
divergences, the counterterms have the form 

AC=^- {Z^ - 1) d,<pdy - ^Ag'(j,\ (5) 

where Ag'^ = (^Z'^Zg — 1^ g'^, and Z^ and Zg are power series in g'^ whose coefficients have 
poles in e. At nonzero temperature, the conventional perturbative expansion also generates 
infrared divergences. They can be removed by resumming the higher order diagrams that 
generate a thermal mass of order gT for the scalar particle. This resummation changes 
the perturbative series from an expansion in powers of g^ to an expansion in powers of 

(g) = 9- 

Screened perturbation theory, which was introduced by Karsch, Patkos and Pe- 
treczky [|14| , is simply a reorganization of the perturbation series for thermal field theory. It 



can be made more systematic by using a framework called "optimized perturbation theory" 



that Chiku and Hatsuda |15] have applied to a spontaneously broken scalar field theory. 



The lagrangian density is written as 

CsPT = -£o + \d,(^d^(^ - ^(m^ - m?)02 _ 1^2^^ ^ ^ ^^^^^ ^ (g) 

where £q is a vacuum energy density parameter and we have added and subtracted mass 
terms. If we set £q = Q and m\ = m?, we recover the original lagrangian @). Screened 
perturbation theory is defined by taking m? to be of order g^ and m\ to be of order g"^, 
expanding systematically in powers of g'^, and setting m\ = m? at the end of the calculation. 
This defines a reorganization of perturbation theory in which the expansion is around the 
free field theory defined by 

£f,ee = -So + \d^<i>d''4> - \rn^'' . (7) 

The interaction term is 
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int 



2/4,1 2 ,2 

q (p H — mi 

24 2 



SPT 



(8) 



At each order in screened perturbation theory, the effects of the term in (H) are included 



to all orders. However, when we set mf 



m 



the dependence on m is systematically 



subtracted out at higher orders in perturbation theory by the ml term in (|^). At nonzero 
temperature, screened perturbation theory does not generate any infrared divergences, be- 
cause the mass parameter in the free lagrangian (|^) provides an infrared cutoff. The 
resulting perturbative expansion is therefore a power series in and mf = whose coef- 
ficients depend on the mass parameter m. 

This reorganization of perturbation theory generates new ultraviolet divergences, but 
they can be canceled by the additional counterterms in A£gpT. The renormalizability of the 
lagrangian in guarantees that the only counterterms required are proportional to 1, 0^, 
df^(j)d^(j), and 0^. With dimensional regularization and minimal subtraction, the coefficients 
of these operators are polynomials in a = g'^/lGn'^ and — The extra counterterms 
required to remove the additional ultraviolet divergences are 



A£gpT = -ASq - - (Am^ - Ami 



(9) 



The vacuum energy counterterm has the form AEq = Ze {m — mi) , where Ze is a power 
series in a whose coefficients have poles in e. The mass counterterms have the form Am^ = 
{Zfj,Zm — ^)m'^ and Am\ = {Z^Zm — ^)ni\, where is the same wavefunction renormalization 
constant that appears in (|^) and Zm is also a power series in a whose coefficients have poles 
in e. 

Several terms in the power series expansions of the counterterms are known from previous 
calculations at zero temperature. The counterterms Ag'^ and Am^ are known to order . 
We will need the coupling constant counterterm only to leading order in a: 



Ag' 



— a + ... 

2e 
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(10) 



We need the mass counterterms Am^ and Am\ to next-to-leading order and leading order 
in a, respectively: 



Am^ 
Ami 



1 

— a ■ 
2e 

1 

— a + ... 

2e 



1 

2e2 24e 

2 



m 



m 



(11) 
(12) 



The counterterm for ASq has been calculated to order [jT6[. We will need its expansion 
only to second order in a and ml: 



{AnfASo 



1 1 

4^ ^ 8e 
1 

— + 



■a + 



48e3 



5 1 
72^ ^ 96e 



a 



m 



- 2 



4e 8e' 



m^m H mi 

^ 4e ^ 



(13) 
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III. MASS PRESCRIPTIONS 



The mass parameter m in screened perturbation theory is completely arbitrary. To com- 
plete a calculation in screened perturbation theory, it is necessary to specify m as a function 
of g and T. One of the complications from the ultraviolet divergences is that the parame- 
ters So, m^, g"^, and mf all become running parameters that depend on a renormalization 
scale II. In our prescription for recovering the original theory, we must therefore specify the 
renormalization scale /i* at which the lagrangian (^) reduces to (0). The prescription can 
be written 



Soiji*) = , 



(14) 
(15) 



where m^{T) is some prescribed function of the temperature. This is the only point where 
temperature enters into SPT. We proceed to discuss the possible prescriptions for m^.{T). 

The prescription of Karsch, Patkos, and Petreczky for m^{T) is the solution to the 
one-loop gap equation: 



-a(/i*) 



Ji{mjT)T^ 



21og^ + l 



(16) 



where the function Ji{x) is defined in ( [A. 8 
weak-coupling limit, the solution to ( 



Their choice for the scale was jj,^ = T. In the 
is m^, = g{^^)T / \/24:. 

There are many possibilities for generalizing (|I6|) to higher orders in g. One class of 
possibilities is to identify with some physical mass in the system. The simplest choice is 
the screening mass iris defined by the location of the pole in the static propagator: 



p2 + + n(0, p) = at p^ 



(17) 



where 11 (po, p) is the self-energy function. Another choice is the rest mass of the quasi- 
particle: rUq = Re a;(0), where uj{p) is the quasiparticle dispersion relation which satisfies 
—u^ + n(i(a; -|- ie),p) = 0. The quasiparticle mass is more difficult to calculate than 
the screening mass. 

Another mass prescription that generalizes (0) to higher orders is to identify m^, with 



the tadpole mass defined by = g {(f) ). This can also be expressed as a derivative of the 
free energy: 



d 



m. 



dm' 



:^{T,g,m,mi,^) 



(18) 



mi=m 



where the partial derivative is taken before setting mi = m. An advantage of the tadpole 
mass IS that (02) is easier to calculate at higher orders than the self-energy 11. 

There is another class of prescriptions that is variational in spirit. The results of SPT 
would be independent of m if they were calculated to all orders. This suggests choosing m 
to minimize the dependence of some physical quantity on m. Taking that physical quantity 
to be the free energy, the prescription is 
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^^•^(7", 9ifJ'),m, mi = m, /i) = . (19) 

We will refer to the solution to this equation as the variational mass. 

One mass prescription that may seem appealing is to choose m* (T) so that the perturba- 
tive approximation is thermodynamically consistent |T^. Given a diagrammatic expansion 
for JF, the entropy density S has a diagrammatic expansion given by 

d 

Sdiag = --^J^iT, g, m, mi, fi) , (20) 

where the partial derivative d/dT is taken with all the other variables g, m, mi, and /i held 
fixed. The entropy density can also be defined by the thermodynamic relation 

^thermo = -^J^{T, g{^i),m = , mi = , fj.) . (21) 

The total derivative takes into account the explicit dependence on T, the T-dependence of 
m^.{T), and also the T-dependence of the running coupling constant if we choose a scale ix 
that depends on T. If the thermodynamic expansions for JF and S were known to all orders, 
there would be no dependence on m or /x, and and ( PT| ) would be equivalent. If the 



diagrammatic expansion is truncated and if any of the parameters g, m, mi, and /i is allowed 
to depend on T, then S may not satisfy (pi]) . An approximation is called thermodynamically 
consistent if S satisfies (^) exactly. This requires 

do dT du dJ^ dm dT dmi dT „ 

I L_ I I i = Q (22) 

dT dg dT dfi dT dm dT dmi ' ^ ^ 

If JF were known to all orders, it would be independent of m and mi at m = mi. Thermo- 
dynamic consistency could then be guaranteed by taking the scale fi to be any function of 
T and choosing g{fi) to be the running coupling constant at that scale. If we only have a 



perturbative approximation to JF, (|22|) is satisfied only up to higher order corrections. One 
way to guarantee thermodynamic consistency is to choose /i = am with a a constant and 
impose the condition 

JF(T, g{am),m, mi = m, (x = am) = . (23) 



dm"^ 



This differs from the variational gap equation (|T9[) only in that we have set fi = am before 
differentiating. This equation does not reduce to the one-loop gap equation ([T6|) at leading 
order, so we will not consider it any further. We will be satisfied by approximations that 
are thermodynamically consistent only up to higher orders in perturbation theory. 

IV. FREE ENERGY TO THREE LOOPS 

In this section, we calculate the pressure and entropy density to three loops in screened 
perturbation theory. The diagrams for the free energy that are included at this order are 
those shown in Fig. ^ together with diagrams involving counterterms. 
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Oa 1a lb 2a 2b 2c 2d 



FIG. 3. Diagrams for the one-loop (Oa), two-loop (la and lb), and three-loop (2a, 2b, 2c and 
2d) free energy. 

A. One-loop free energy 



The free energy at leading order in g"^ is 



(24) 



where Aq^^o is the term of order in the vacuum energy counterterm (p!3[) . The expression 
for diagram Oa in Fig. H is 



Oa 



1 

2jp 



log + 



(25) 



The sum-integral in (^) is over the Euclidean momentum P = (ct;„, p) and we define 
P^ = The sum-integral includes a sum over Matsubara frequencies Un = 27mT and 

a dimensionally regularized integral over the momentum p with a measure that is defined in 
Appendix |^. In dimensional regularization with 3 — 2e spatial dimensions, the diagrams for 
JF have dimensions (energy)^" To obtain the renormalized free energy density with dimen- 
sions (energy)"^, we multiply the diagrams by /i^^, where fi is an arbitrary renormalization 
scale, before taking the limit e — 0. The coupling constant in dimensional regularization 
is gfi"^, where g is the dimensionless renormalized coupling constant. Including the overall 
factor of /i^^ and the factor of /i*^ from the coupling constants, there is a factor of /i^'^ for 
each sum-integral. We choose to absorb this factor into the measure of the sum-integral. 

The sum-integral in ( pSf ) is expressed as a function of e in the Appendix. It has a pole 
at e = 0. The result for the diagram is 



Oa 



4(47r) 



fx 



m 



2e 



1 3 21 + 7r2 45 + 37r2 + 4< 1 , 

- H \ e H Z_L^g^ 

e 2 12 24 



(26) 



where Jq is the function oi m/T defined in (|A.5|) . We have kept all terms that contribute 
through order e^, because they enter into higher order diagrams involving counterterms. The 
pole in e in (|26| ) is canceled by the zeroth order term Aq^^o in the counterterm (p!3|). The 
final result for the one-loop free energy is 



(47r)2^o = (4vr)2^o - ^(2^ + 3)^^ - \JoT' , 



(27) 



where L = log(/i^/m^) and Jo can now be replaced by its value at e = 0, which is given 
in (KM- 
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B. Two- loop free energy 



The contribution to the free energy of order is 



= ^la + ^ib + ^i£o + ^-^Aim^ 



(28) 



where AiSq and Aim^ are the terms of order g"^ in the counterterms ( pT]) and ([T3|), respec- 
tively. The expressions for the diagrams la and lb in Fig. |^ are 



^ la 



1 



p P"^ + m? 



1 



2 "~'Tp P'^ + m? 
The results for the diagrams can be expressed as 



(29) 
(30) 



^ la 



^(47r)2 \m) 



1 2 18 + 7r2 + + ' 
+ ^ e 



e 

-2 



6 



m 



1 12 + 7r2 ■ 

e 12 



lb 



2(47r)2 Vm 



2e 



1 , 12 + 772 

- + IH e 

e 12 



+ Jit4 , 



(31) 
(32) 



where a = g'^/lGir'^. We have kept all terms that contribute through order e, because they 
are needed for counterterm diagrams in the three-loop free energy. The poles in e in (|3lD 
and ( |32D are canceled by the counterterms in (|28|). The final result for the two-loop free 
energy is 



{Anfj^i = -\{L + l)m2 - JiT^l + -a \{L + l)m'^ - J^T^'^ 
2 L -I 8 '- 



(33) 



C. Three-loop free energy 

The contribution to the free energy of order g^ is 

^2 = ^2a + ^2b + ^2c + ^2d + A^So + ^^^m^ + {A.m'f 

y om^ om^ J g^ mf 

where we have included all the appropriate counterterms. The expressions for the diagrams 
2a, 2b, 2c, and 2d in Fig. 13 are 
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^ 2a 
^2b 



-16^ 



-.9 



48'' TPQR (P2 + rn2)(g2 + m2)(^2 + ^2)((p + g + /2)2 + ^2 



^2c 4^?X^ P^ + m^tq (Q2 + ^2)2 



^2d = -^ml 



1 



4 (P2 + m2)2 • 

The results for these diagrams in the hmit e ^ are 



a 



2a 



16(47r)2 \mj 



6e 



1 ^ 2 ^ 12 + ^ 8 + vr^ + ^/>"(1)' 



+ 



4e 

1 2 18 + 7r2 
- + - + 



m 



6 



Jam^ - 2 



1 1 6 + 7r2 
- + - + 



6 



2r7-i2 



1 

-2 - + 1 

e 



2b 



'48(47r)2 Vm/ 



6e 



2 23 35 + 7r2 ^' 
+ — + — ^ + Co 



6 



3e2 



12 



2e 



4 

m — 



6 17 

- + 4Ci 



2t^2 



Jim^T 



JfT^ + [QK2 + AKs] 



am 



2c 



4(47r)2 Vm) 



4e 



1 1 6 + 

-T + -^ 



6 



ml 



2d 



2e 



fl 



4(47r)2 Km) I e ^ ' ' 



(35) 
(36) 
(37) 
(38) 



(39) 

(40) 
(41) 

(42) 



The poles in e are canceled by the counterterms in (|34D . The final result for the free energy 
is 



(47r)2j-2 = -1 (L + J2) - ^(L + J2) [(L + l)m^ - J,T^ 



mi 



——a 
48 



(SL^ + ITL^ + f L - 23 - ffTT^ - ^"(1) + Co + 3(L + 1)^2) m' 
- (12L2 + 28L - 12 - TT^ - 4Ci + 6(L + 1)72) Jim^T^ 
3(3L + 4) J2 + 3J1V2 + QK2 + 47^3) T 



(43) 



D. Pressure to three loops 



The pressure V is given by —!F. The contributions to the free energy of zeroth, first, and 
second order in g"^ are given in (|27|) , (|33|), and (|^), respectively. Adding them and setting 



£^0 = and ml = m^, we get the approximations to the pressure in screened perturbation 
theory. The one-loop approximation is obtained by setting £^0 = in 
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(47r)2po = - 4 JoT^ + (2L + 3) 



(44) 



The two-loop approximation is obtained by adding (|33D with ml = m^: 



(47r)2po+i 



4JoT^ + AJim^T^ - {2L + 1) m 



-a 



The three-loop approximation is obtained by adding ( ^31 ) with m\ = m?: 
1 



(47r)2po+i+2 



4 JqT* + AJim'T^ + 2 Jam^ - 



— a 
8 

H 

48 



(45) 



JiT^ - (L + l)m2] [jiT^ + 2J2m2 + (L - l)m2' 

3J2 (jiT^ -{L + + (3(3L + A)Jl + 67^2 + AK^) 

- (l2L^ + 28L - 12 - TT^ - 4Ci) Jim^^T^ 



+ (5L3 + ITL^ + f L - 23 - ffyr^ - ^"{l) + Co 



(46) 



where L = log(/i^/m^), Cq = 39.429, Ci = —9.8424, the J„'s are the functions oim/T given 
in ([A.8| ), and K2 and are functions of m/T given in Ref. |]18[. Note that the dependence 
on L has canceled from the term proportional to a° in (^). 



E. Entropy to Three Loops 



The perturbative expansion for the entropy density S is defined in (|2l|) . The one-, two-, 
and three-loop approximations to S are obtained by taking the partial derivatives with re- 
spect to T, with a, m, and /x fixed, of the expressions for the pressure in (^4]), (|45|), and (^61). 
The partial derivatives of the functions J„(/3m) can be evaluated using the recursion rela- 
tion ( |A.6| ). The partial derivatives of Kn{(3m) can be evaluated numerically. 



The one-loop approximation is obtained by differentiating (|4^) : 



The two-loop approximation is obtained by differentiating (^5]) 



'0+1 



2 JoT* + 2 Jim^T^ + Jgm' 



1 

a 

2 



JiT^ -{L^ l)m^ JiT^ + Jam 



The three-loop approximation is obtained by differentiating (^61): 
1 



(47r)2T5o+i+2 



1 



2 



+J3 f^iT' - (i^ + l)m') 



(47) 



(48) 



12 



+ 24" 



+6J2 (jiT^ - (L + l)m2) (jiT^ + J^m 



+ (6(3L + 4) J2 + 127^2 + SfsTs) - (37^2 + 2/^3) mT 
+6{3L + A)JiJ2m'^T^ 

- (l2L^ + 28L - 12 - TT^ - 4Ci) (jiT^ + J^m 
The primes on and denote differentiation with respect to (3m. 



m 



(49) 



V. SCREENING MASS TO TWO LOOPS 

In this section, we calculate the screening mass to two loops. The diagrams for the self 
energy that are included at this order are those shown in Fig. ^ together with diagrams 
involving counterterms. The screening mass is the solution to the equation ([171) . This 
equation can be solved order- by-order in powers of a and mf. The solution at zeroth order 
in g"^ is simply = m^. 




1a 2a 2b 2c 



FIG. 4. Diagrams for the one-loop (la and lb) and two-loop (2a, 2b and 2c) self-energy. 



A. One-loop self-energy 



The self-energy at first order in g'^ is 

Hi = Hia -ml + Aim^ 



(50) 



where Aim^ is the mass counterterm of order a given in (pTf) . The expression for the diagram 
la in Fig. H is 



Hie 



P2 + 



m 



2 ■ 



The result for the diagram is 



(m) 


'_ 







1 , 12 + 7r2 

- + IH e 

e 12 



(51) 



We have kept all terms that contribute to order e, because they are needed for counterterm 
diagrams in the two-loop self-energy. The pole in e in (|5TD is canceled by the counterterm 
Aim^. The final result for the one-loop self-energy is 
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(52) 



B. Two- loop self- energy 



The contribution to the self-energy of second order in g"^ is 



MP) = M + MiP) + + f^Aim^ + ^A^g^ + A^m^ - A.mj . 



9^ 



(53) 



The expressions for the diagrams 2a and 2b in Fig. 11 are 



= -Tfl- 



n2b(P) 



4^ Tq Q"^ + m?TR {R^ + m?f ' 
6^ Tqr Q'^ + m?m + m?{P + Q + RY + w? 



(54) 
(55) 
(56) 



The diagrams and n2c are independent of the momentum P. The resuhs for these 
diagrams in the hmit e — >^ are 



4 \m, 

n2c = — 

2 Vm 



' ' ' ' ^-±^)m^+(l + l)j2m^-(l + J2)j.T^ 



e2 e 



2e 



+ J2 



(57) 
(58) 



The diagram n2b depends on the external momentum P. The equation ([TTD for the screening 
mass involves the self-energy at po = 0. To calculate the screening mass to second order 
in g"^, we need the analytic continuation of 11(0, p) to = — m^. This is calculated in the 
appendix. The result is 



n2b(o,p) 



a 

p^=— m-^ 



4e 



3 17 



m — 3 



rp2 



(59) 



The poles in (p7D-(^) are canceled by the counterterms in (|55D. The final result for the 
two-loop self-energy at po = and p^ = — is 



n2(o,p^ 



1 1 1 r 

-a{L + J2)ml + :t7«^S 12L2 + 28L - 12 - vr^ - 4Ci + 6 (L + 1) J2 

^ I 



m 



(3L + J2) Ji + 2i^i + 2ir2 



(60) 
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C. Screening mass 



Since the dependence of the self-energy on the momentum enters only at order g'^ and 
since the leading-order solution to the screening mass is m<j = m, the solution to the equa- 
tion (0) to order g'^ is simply 



= + n(0, 



(61) 



We proceed to calculate the expression to order g'^ and to order g'^. 

The solution to order g^ is obtained by inserting the one-loop self-energy ( p2D into (^TJ). 
Setting ml = m?, the result is 



= -a 



JiT^ -(L + l)m2 



(62) 



If we choose m = mg = m^, this is identical to the one- loop gap equation ([T6|). 

The solution to order g^ is obtained by inserting the sum of ( ^21) and (|60|) into (^ 
Setting m\ = m^, the result is 



= -a 
' 2 



JiT" + (J2 - iW 
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-a 



6J2 ( JiT^ -iL + l)m^ 



+6 (sLJi + 2ki + 2K2) - [l2L^ + 28L - 12 - vr^ - 4Ci 



m 



(63) 



Note that the dependence on L has canceled in the order-a terms. 



VI. GAP EQUATIONS 

In this section, we solve the gap equations that determine the arbitrary mass param- 
eter in screened perturbation theory. We consider the one-loop gap equation and three 
generalizations to a two-loop gap equation. 



A. One-loop Gap Equation 



The one-loop gap equation is given in (|T6|) . It is convenient to introduce the gap function 
defined by 



G = - -« [jiT^ - (L + l)m^ 



(64) 



The one-loop gap equation is then G = 0. For simplicity of notation, we will often suppress 
the subscripts * on m and fi. 

Before solving the one-loop gap equation, we need to choose a value for fi. It is natural 
to take yU to be proportional to one of the two energy scales in the equation, T and m. We 
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will consider two possibilities, /i = a(27rT) and /i = am, and allow the coefficient a to vary 
from I to 2. Given either of these choices for /x, the gap equation can be solved for m as a 
function of Q;(/i). The renormalization group equation can then be used to express a(/i) 
as a function of aiT.nT). 

In the weak-coupling limit fif — > 0, the solution to the gap equation G = approaches 



27r^ 



a(/i)T^ 



1 - y6a^/2 - ( lo; 



+ 7 - 3 I a + 0{a 



3/2> 



A-nT 

In the strong-coupling limit g ^ oo, the gap equation reduces to 

.2 1 



21og^ + l = 
m 

This has a solution only if /x > e 



X 



-1/2 







m. 



(65) 



(66) 



E 0.7 



One Loop 

nJ <[i< 4nT 
■ ^nn*<n<2m^ 




FIG. 5. Solutions m^,{T) to the one- loop gap equation (shaded bands) and the two- loop gap 
equations (lines) as functions of g{2-KT). 

In Fig. ^, the solutions m^:{T) to the one-loop gap equation as a function of g{2iTT) are 
shown as bands obtained by varying /i by a factor of two around the central values = 2ttT 
and /i = m*, respectively. The solutions are normalized to the leading-order screening mass 
rriLo = g{2TTT)T/V2A. 



B. Screening Gap Equation 

The screening gap equation is obtained by identifying m with m^. The one- loop expres- 
sion for the screening mass is given in (0). Thus the one-loop screening gap equation is 
simply G = 0. The two-loop expression for the screening mass is given in (^). The two- loop 
screening gap equation can be written as 
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1 - 2« (-^2 + L) 



12 (l Ji + + k^) 



6L^ + 22L - 12 - TT^ - 4Ci m' 



. 



(67) 



From this expression, it is easy to see that the solution m to the gap equation differs from 
the solution ( |65| ) to the one- loop gap equation by terms of order a^T^. The weak-coupling 
expansion of the solution must of course agree through order a^T^ with the weak-coupling 
expansion of given in (|^). 

The solutions to the screening gap equation for fi = 2ttT and /i = are shown in 
Fig. |. In the case /x = 27rT, the screening gap equation cannot be continued beyond 
g{2'KT) = 2.60. For fi = irT, it terminates at g{27TT) = 2.31, while for /i = 4nT, it 
terminates at g{2nT) = 3.04. If we choose /i = am^,, the solution can be continued to much 



larger values of g. For /i 
with fi = m^,. 



it lies very close to the solution to the one-loop gap equation 



C. Tadpole Gap Equation 



The tadpole mass is defined in (|T^). The one- loop expression is given by differen- 



tiating ([27|). The result is identical to the one-loop expression (|62D for the screening mass. 
To obtain the two-loop expression for the tadpole mass, we add the one- and two-loop free 
energies (pTj) and (p3D, differentiate with respect to m^, and then set m\ = w?. The result 
is 



JiT^ + {J2 - iW] - -a" (J2 + L) [jiT^ -{L + l)m^ 



(68) 



The order-a term is identical to that of the screening mass (|63|), but the order-a^ term is 
much simpler. 

The one-loop tadpole equation is simply G = 0. The two-loop tadpole gap equation is 
obtained by setting rrit = m in ([681). It can be written in the form 



1 - -a ( J2 + L) 



G = . 



(69) 



Thus the two-loop tadpole gap equation is identical to the one-loop gap equation: G = 0. 
The solutions for /x = 2ttT and n = m^, are at the centers of the shaded bands in Fig. ^ 



D. Variational Gap Equation 



The variational mass m„ is the solution to (|T^). The one-loop variational gap equation 
is obtained by differentiating the two-loop expression ( ^Sj) for the pressure with respect to 

and setting it equal to zero. This gives (L -|- J2) rn?G = 0, which reduces to the one-loop 
gap equation: G = 0. 

The two-loop variational gap equation is obtained by differentiating the three-loop ex- 
pression (H^) for the pressure. It can be expressed in the form 
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\a{J, + LfG-\U 



+ 48" 




- 12{L + 2) Ji J2 + 



f5m 



(70) 



where K2 and K'^ are the derivatives of K2 and with respect to (3m. In the coefficient of 
a^T^, we have written exphcitly only the terms that are singular as (3m — > 0. The l/{(3m)^ 
singularities cancel between the Jf and K2 term. If we keep the most singular terms in the 



coefficients of each of the three terms in (|70D , the equation reduces to 



■G 



TT 



(/3m) 2 4 (/5m) 3 
- 32n%L + 2) - (3A;2 + 2^:3) log(/?m) 



3(A;2 + A;^) - 2{h + 4) 



2m2 



48(3m 







(71) 



where /c2log(/?m) + k2 and A;3log(/3m) + are the coefficients of (3m in the small-/3m 
expansions of K2 and i^^a, which are given in ( |A.16| ) and ( [A.ITI ). 

The solution to the quadratic equation ([7T|) for G is proportional to a(3mT'^. The 
solution to the gap equation therefore differs from the solution ( |65D to the one-loop gap 
equation by terms of order a^/^T^. This is a little disturbing, but even more disturbing is 
the fact that ( |7TD has no real-valued solutions for G unless L < 2.0984 log(/?m) + 4.1541. If 
we assume that m gT /\/2A as (7 — 0, then this condition is violated for sufficiently small 
g whether we set /x = a(2'7rT) or /i = am. Since there are no solutions in the neighborhood 
of (yf = 0, we will not consider the two- loop variational gap equation any further. 



VII. SPT-IMPROVED OBSERVABLES 

In this section, we use the solutions to the gap equation in Sec. ^ to obtain successive 
approximations to the pressure, screening mass, and entropy in screened perturbation theory. 

A. Pressure 



The two-loop SPT-improved approximation to the pressure is obtained by inserting 
the solution to the one- loop gap equation (|16D into the two- loop pressure (^). We can 
simplify the expression by using (p!6| ) to eliminate the explicit appearance of logarithms of 
/i. Remarkably, this eliminates all the terms of order a and the expression reduces simply 
to 



(4vr)2p, 



0+1 



4JoT^ + 2 Jim^T^ + m 



(72) 



The Jo term in (|72D is the pressure of an ideal gas of particles of mass m. Inserting the 
solution to the one-loop gap equation shown in Fig. ^, we obtain the bands shown in Fig. |^. 
The lower and upper bands correspond to varying yU by a factor of 2 around the central 
values /i = 27rT and /i = m^,, respectively. 
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— - l^ = m* 



0.94 I ■ ' ■ ' ■ ' ■ 1 

12 3 4 

g(27rT) 

FIG. 6. Two-loop (shaded bands) and three-loop (lines) SPT-improved pressure as a function 
of 5(27rT). 

The three-loop SPT-improved approximation to the pressure is obtained by inserting 
the solution to a two- loop gap equation into the three- loop pressure (^61) . In Fig. we 
show the three- loop SPT-improved pressure as a function of g{2TTT) for different two- loop 
gap equations. The solid hne is the result using the two-loop screening gap equation with 
/i = 27rT. It cannot be extended past g{27TT) = 2.60. The dashed line is the result using 
the two-loop tadpole (or one-loop) gap equation with /i = 27iT. The dotted line is the result 
using either the two-loop screening gap equation with /x^, = or the two-loop tadpole 
gap equation with /i^, = m^,. The two are indistinguishable on the scale of the figure. The 
variations among the three-loop SPT-improved approximations for the pressure are much 
smaller than one might have expected from the variations among the screening masses. For 
example, at g{2nT) = 2, the solutions to the two-loop gap equations shown in Fig. |^ vary by 
about 12%, while the three-loop approximations to the pressure shown in Fig. ^ vary only 
by about 0.07%. 

Since the solution to the screening gap equation at = a(27rT) cannot be continued 
beyond a critical value of g and the solution for fi = am^, is close to the solution to the 
tadpole gap equation for /i = am^,, we will consider only the tadpole gap equation from now 
on. In Fig. |^, we show the one-, two-, and three- loop SPT-improved approximations to the 
pressure using the tadpole gap equation. The bands are obtained by varying /z by a factor 
of two around the central values = 2ttT and /i = m^. The one- loop bands in Fig. |^ lie 
below the other bands; however, the two- and three-loop bands all lie within the g^ band 
of the weak-coupling expansion in Fig. |l]. The one-, two-, and three-loop approximations to 
the pressure are perturbatively correct up to order g^, g^, and g^, respectively; however, we 
see a dramatic improvement in the apparent convergence compared to the weak-coupling 
expansion. 

The choice /i = am^, appears to give better convergence than /i = a(27rT), with the 
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(a) TiT < n < 47iT 



(b) im*<^< 2m, 
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FIG. 7. One-, two-, and three-loop SPT-improved pressure as a function of g{2TTT) for (a) 
ttT < /j, < AttT and (b) < ii < 2m*. 



three-loop band falling within the two-loop band. The bands for jj = am^ are narrower 
than those for yU = a(27rT) partly because /i = a{27iT) is larger and therefore closer to the 
Landau pole of the running coupling constant. If g{2TTT) = 2, the Landau pole associated 
with the five- loop beta function is far away at fi = 2.11 x 10^(27rT). If g{27TT) = 4, the 
Landau pole is rather nearby at /i = 5.49(27rT). The coupling constant g{m^) is smaller 
than g{2TTT), having the values 1.76 and 3.07 if g{2TTT) = 2 and 4, respectively. Choosing 
II = am^: instead of /i = a(2'7rT) will therefore make the error due to the terms in the 
pressure smaller by factors of about 0.60 and 0.35 respectively. The band m^../2 < fi < 2m^ 
may therefore give an underestimate of the error of SPT. 



B. Screening Mass 



The one-loop SPT-improved approximation to the screening mass is simply the 
solution m^{T) to the tadpole gap equation. A two- loop SPT-improved approximation can 
be obtained by inserting the solution to the gap equation for m into (|63D . In Fig. ||, we show 
the one-loop and two-loop SPT-improved approximations to the screening mass as functions 
of g{2iTT). The bands are obtained by varying /i by a factor of two around the central values 
fi = 2ttT and /i = m,,. 

The choice fi = am^ appears again to give better convergence than fi = a{2TrT), with 
the two- loop band falling within the one- loop band. With fi = am^, there is a dramatic 
improvement in apparent convergence over the weak-coupling approximations, which are 
plotted on the same scale in Fig. However, there is not much improvement in the apparent 
convergence with fi = a{2nT). The conservative conclusion is that screened perturbation 
theory is not as effective in improving the prediction for the screening mass as it is for the 
pressure. 
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(a) nJ<[i< 4nT 
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FIG. 8. One-loop and two-loop SPT-improved screening mass as a function of g(2TTT) for (a) 
ttT < h < AttT and (b) ^m^, < ^ < 2m*. 

C. Entropy 

The one-, two- and three-loop SPT-improved entropies are obtained by replacing m in 
the expressions (^)-p9|) for 5o, i5o+i, and iSo+1+2 with the solution m^, to the one-loop gap 
equation G = 0. Using the gap equation to eliminate the logarithm L, the expression for 
the two-loop entropy reduces to 



(47r)"T5o+i = 2JoT^ + Jim^T^ 



(73) 



This is identical to the one- loop expression (0), which is the entropy of an ideal gas of 
particles with mass m. In Fig. ^, we show the two- and three-loop SPT-improved approxi- 
mations to the entropy as functions of g{2TTT). The entropy density is normalized to that of 
an ideal gas: iSideai = (27r^/45)T^. The bands in Fig. |^ correspond to varying /i by a factor 
of two around the central values 27rT and m^,. Once again, the choice /i = am^, seems to 
give better convergence with the three-loop band lying very close to the two-loop band. 

The entropies shown in Fig. ^ are successive approximations to the diagrammatic en- 
tropy defined by Eq. However, the entropy can also be defined by the thermodynamic 
relation (^If ). Thus successive approximations to S can be obtained by differentiating the 
pressures shown in Fig. ^ with respect to T. In that figure we show the ratio of the pressure 
to that of an ideal gas as a function of g{27TT). Defining the function f{g) by 



V{T) = Pidcai(T)/(5(27rr)) , 



the thermodynamic entropy is then given by 



<Sthermo(T) — Sideal(T) 



fig) 



27r2 



(74) 



(75) 



where g = g{2'7iT), a = g^/lGn'^, and /?(a) is the beta function given by the right side of 
(§). In Fig. |10|, the black curves are the two- and three- loop diagrammatic entropies for 
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FIG. 9. Two-loop and three-loop SPT-improved entropy as a function of g{2TTT) for (a) 
ttT < h < 4ttT, (b) < ^ < 2m*. 

/i = 27rT and yU = m^:. The gray curves are the the thermodynamic entropies obtained from 
the one-, two-, and three-loop SPT-improved pressures. One can see clearly the approach to 
thermodynamic consistency as one goes from the two-loop to the three-loop approximation. 
With the choice fi = 27rT, the two-loop entropy is almost perfectly thermodynamically con- 
sistent. However, this is probably an accident because the deviations from thermodynamic 
entropy are evident in the three- loop entropy. With fi = am^,, deviations from thermody- 
namic consistency are very small for both the two- and three-loop entropies. This is another 
indication that SPT improvement is more effective if we take the scale /i to be much lower 
than 27rT. 
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FIG. 10. SPT-improved entropy as a function of g{2-iTT) compared to the thermodynamic 
entropy for obtained from the SPT-improved pressure for (a) /i = 2ttT and (b) fi = m*. 
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VIII. CONCLUSIONS/DISCUSSION 



We have studied the effectiveness of screened perturbation theory in reorganizing the 
perturbation series for a thermal scalar field theory. We applied it to the pressure and the 
entropy calculated to three loops and to the screening mass calculated to two loops. 

We considered three alternatives for generalizing the one-loop gap equation to two-loop 
order. The most useful turned out to be the tadpole gap equation, which at two loops is 
identical to the one-loop gap equation proposed by Karsch, Patkos and Petreczky. The 
solution to the two-loop variational gap equation does not match onto the one-loop gap 
equation in the weak-coupling limit. The solution to the two-loop screening gap equation 
cannot be extended above g{2TTT) = 2.60 if we choose the scale to be /i = 2ttT. 

The predictions of SPT depend on an arbitrary scale fi that arises both from the renor- 
malization of the coupling constant and from the renormalization of ultraviolet divergences 
introduced by screened perturbation theory. These two effects could be separated by in- 
troducing two renormalization scales, /is and /X4. These scales would be associated with 



contributions from soft and hard modes respectively as in Ref. |]19[. One way disentangle 
the dependence on these scales would be to evaluate the integrals as expansions in m/T. We 
evaluated our integrals by integrating numerically over all momenta, which precluded any 
separation of the scales. Instead, we considered two possibilities for the scale, ^ = and 
/i = 27rT, which correspond to the central values expected for /is and /i4, respectively. We 
allowed for variations of /i around these central values by factors of two to provide a lower 
bound on the theoretical uncertainty. The choice /i = gives smaller bands from varying 
the scale, but this is largely due to the fact that the coupling constant gim^) is smaller 
than g{2TTT). Thus the size of the bands is not a good indicator of the success of the SPT 
improvement. 

A better indication of the success of SPT improvement is the stability of the predictions 
as you go to higher order in the loop expansion. The choice = 2ttT gives a significant 
improvement in stability for the pressure compared to the weak-coupling expansion. However 
the SPT improvement seems to be much more effective using /i = m^, than /i = 2tiT . The 
three-loop band lies within the two-loop band for the pressure and it lies very close for the 
entropy. The two-loop band also lies within the one-loop band for the screening mass. The 
two-loop and three-loop approximations for the entropy are also very close to thermodynamic 
consistency if we choose /i = m^,. If we set ^ = 2ttT, then going from the two-loop to the 
three-loop approximations to the pressure or entropy moves the prediction closer to that for 
II = m^. This suggests that the SPT-improved prediction for /x = m^, is more accurate than 
that for /i = 2ttT. All this evidence indicates that SPT improvement is most successful if ix 
is taken to be much smaller than 2'kT . 

To remove the additional ultraviolet divergences introduced by SPT, we have chosen 
to use dimensional regularization with modified minimal subtraction. This choice is of 
course not unique. For example, we could have also chosen to subtract the piece of the free 
energy that is independent of T for fixed m, i.e. Tb. = ^{T, g,m, fx) — T{T = 0,g,m,fi). 
This would result in a different reorganization of the perturbation series that would also 
agree with the exact result if summed to all orders. We take into account the theoretical 
uncertainty associated with the choice of subtraction scheme by allowing variation of the 
renormalization scale fi. For example, by examining Eqs. (|4^) and (^), we can see that 



23 



the alternative described above corresponds to L = —3/2 or the one-loop approximation 
and to some value in the range — | < L < — 1 for the two- loop approximation. Setting 
/i = am* with | < a < 2 corresponds to varying L in the range —1.4 < L < 1.4. Therefore, 
this variation of /i does take into account the ambiguity associated with the subtraction 
scheme at T = It would of course be preferable to separate the ambiguity from the SPT 
subtractions from the ambiguity from renormalization of the original theory by allowing 
separate renormalization scales /is and fi^ as mentioned above. 

Our results demonstrate the effectiveness of screened perturbation theory in providing 
stable and apparently converging predictions for the thermodynamic functions of a massless 
scalar field theory. An essential ingredient of this approach is using the solution to a gap 
equation as the prescription for the mass parameter m. This success of screened perturbation 
theory adds support to the proposal of Ref. |]19| to use HTL perturbation theory to reorganize 



the weak-coupling expansions for the thermodynamic functions of QCD. In Ref. |]T9|, the free 
energy was computed only to one-loop order, so there was little alternative to using a weak- 
coupling expression for the thermal gluon mass parameter. However, our experience with 
SPT indicates that the stability of the predictions is greatly improved by using a solution 
to a gap equation for the mass parameter. A gap equation can be derived from the free 
energy calculated to two-loop order in HTL perturbation theory. Until that calculation is 
carried out, quantitative comparisons of the predictions of HTL perturbation theory with 
the nonperturbative results of lattice gauge theory are probably premature. 
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APPENDIX A: SUM-INTEGRALS 

In the imaginary-time formalism for thermal field theory, a boson has Euclidean four- 
momentum P = {po,p), with = pI + p^. The Euclidean energy pq has discrete values: 
Po = 27rnT, where n is an integer. Loop diagrams involve sums over po and integrals over 
p. We use dimensional regularization to regularize ultraviolet or infrared divergences. Our 
choice for the measure in the sum-integrals is 

where 3 — 2e is the dimension of space and n is an arbitrary momentum scale. The factor 
{e'^ /Any is introduced so that, after minimal subtraction of the poles in e due to ultraviolet 
divergences, n coincides with the renormalization scale of the MS renormalization scheme. 
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1. One- loop sum-integrals 



The one-loop sum-integrals that appear in the free energy can be separated into a 
temperature-independent term and a term that depends explicitly on T: 
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The thermal terms can be expressed as integrals involving the Bose-Einstein distribution 
function: 

^4-2n-2. 1 
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These integrals satisfy the recursion relation 

xJ'^{x) = 2eJn{x) - 2x^Jn+i{x) . 
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(A.5) 



(A.6) 



The temperature-independent terms in (|A.2|) -( [Al4D can be expanded as a Laurent series 
around e = by using 

e^T(l + e) = l + ^e' + ^^{1^ + Oie') . (A.7) 

The functions J„ have Taylor expansions around e = 0. They often appear multiplied by 
poles in e, but it is counterproductive to expand J„ in powers of e, because the poles always 
cancel in physical quantities. 

If we set e = 0, the integrals for = 0, 1, 2 reduce to 

^4-2n I 



^4-2n 



dk 



lo (p + m2)i/2e/3(fc^W)V2 _i • 
The integral J3 requires a subtraction to remove a linear infrared divergence: 
J3(/5m) = -2/3-2 / dk 

In the limit jSm 



F ^2)1/2 e/3(fe2+m2)i 

0, these integrals reduce to 



/2 



1 mef^"" -1 





45 ' 




Jl - 


47r2 
3 


47r/?m - 


J2 - 


27r 

[3m 




J3- 


TT 


1 


(/3m)3 


(/5m)2 



47r 



+ 7 (/5m)% 



+ 



Atx' 



(A.8) 

(A.9) 

(A. 10) 
(A.ll) 

(A.12) 
(A.13) 
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2. Basketball sum-integral 



The only nontrivial sum-integral required to calculate the free energy to three loops is 
the massive basketball sum-integral: 

^ (A. 14) 



ball 



fpQR (P2 + m2)(Q2 + rn2)(i?2 + ^2)[(p + g + i?)2 + ^2] • 
This sum-integral was evaluated in Ref. pO] and the result is 



ball 



(4vr) 



2 23 35 + 7r2 ^' 
+ — + — + Co 



3e2 



2e 



4 I 



6 17 



+ 4Ci 



6 



+ ( - + 12 ) JtT"" + {QK2 + 47^3) [>T^ , 



(A.15) 



where Cq = 39.429, Ci = —9.8424, and K2 and are functions of /5m. They are expressed 
in Ref. [jl8[ as three-dimensional integrals that can be evaluated numerically. Their behavior 
in the limit f3m — is 

3277-4 

K2 — > —— [log(/3m) - 0.04597] - 372.65 [log(/5m) + 1.4658] /3m , (A. 16) 

9 

K3 — > 453.51 + 1600.0 [log(/3m) + 1.3045] pm . (A.17) 

The leading terms are given analytically in Ref. [|18|. The terms proportional to jSm were 
determined numerically. 



3. Sunset sum-integral 



The only nontrivial sum-integral required to calculate the self-energy to two loops is the 
sunset sum-integral, which depends on the external four-momentum P = {po,p)'- 

1 1 



'fQR Q^ + m^R^ + m^P + Q + Rf + tu^ ' ^^'^^^ 

This sum-integral can be separated into terms with zero, one, and two thermal distributions, 
respectively At po = 0, it can be written as 

Xsun(0, p) = Xi°) (P') + 3XW (P') + 3Xi£(p^) , 



where 



1 



1 



QR Q'^ + m?R^ + m?{P + Q + Rf + m? 



P=(0,p) 



^''Z'"^^^^ lRR^ + m^P + Q + RY + m2 
Xi^l(p') = Re / n5{q) f n5{ ' 

Jn Jr 



(P+Q)2 = -[R2_(p+q)2+ie] 



(p + g + r)2 — m? + ie 



(A. 19) 

(A.20) 
(A.21) 
(A.22) 



p={0,p) 



The integral denotes the dimensionally regularized integral over the Minkowski momentum 

(go, q), and n5{q) = n{qo)2Tr6{q'^ — w?). 
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a. Zero thermal factors 



To calculate the screening mass to two loops, we need the analytic continuation of 



the integrals ( |A.2CI|) - ([A.22|) to the point = —m?. The integral (|A.2CI|) was evaluated in 



Ref. 18 



(4vr) 



4e r 



3 17 
"2^"4^ + ^^ 



where Ci = -9.8424. 



(A.23) 



h. One thermal factor 



The integral ( |A.21|) can expressed as 
1 //i\^^rl 



(47r)4 Vm, 

-8rc^g^^r^x/log 
JO E„ Jo \ 



— x{l — x){E^ — k'^ 



(A.24) 



where k = |p + q| and (...) denotes the angular average. After averaging over angles, ([A.24| ) 
can be analytically continued to = — m^. The result is 



/i 



4e r 



(47r)4 Vm 

where Ki, which is a function of (3m only, is defined by 



Ki = -— / dq 



T2 Jo 



The function fi{x,q) in the integrand is 



' dx fi{x,q) . 



~ x^ + {l-xf ( 2x{l - x)q/m\ 
Ji{x, q) = ^r—r. ^ — ; — atan | — | — 1 



2x(l — x)q/m 



x^ + [1 — x) 



1 , 



+ (1 - xf ' + 4x^(1 - xfc^jm^ 



(A.25) 



(A.26) 



(A.27) 



c. Two thermal factors 



The integral ( |A.22|) can be expressed as 



32 , q^niEg) 

dq 



(47r)4 Jo 



En 



-1] 



E,. 



r (J 



^2 _ _ _|_ 



(A.28) 
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where E^r = Eq + aEj., k = |p + q + r|, and a is summed over ±1. After performing the 
angular average, (|A.28 ) can be analytically continued to = — m^. The result is 



J(2) 



-m 



{Any 



where K2, which is a function of /5m only, is defined by 



''' = T-L E, 



dr 

Ej. 



(A.29) 



(A.30) 



The function in the integrand is 

[E'^ - (g + rff + Am^{q + rf 



f2{E,q,r) = log 



[E'^ — (g — r)^]^ + 4m'^{q — r) 



2(q + r) 2m(q + r) 2\q — r\ 2m|g — r| 

-atan-^^;; ; — : — ttt H atan- 



m 



E^-{q + ry 



m 



E2 - {q-r f 



8mqr\^ E"^ — 



H V E'^ — m^atan — ^ nx / a 

m E4-2(E2-2m2)(g2 + r2) + (g2_r 



2^2 



(A.31) 



If E'^ < m^, the last term in (|A.31|) should be replaced by a manifestly real- valued expression 
using the identity 2ix atan(ia;/|/) — ^ a;log[|j/ — x\/\y + x\]. 

Our final result for the sunset sum-integral evaluated at po = and = 
obtained by combining ( [A:23|) , (^^, and ( [A:29|) as in ( |Al9| ): 



is 



^sun(0,p) 



(47r)4 Vm 



/i 



4e 



17 



+3 



2e2 4e 

1 

e 



(A.32) 



The behavior of the functions Ki and in the limit jSm — > is 



K2 



log 



f3m 



3 + 



C(-i) 
C(-i) 



-47r' 



(A.33) 
(A.34) 



The result (|A.33| ) was computed analytically. The result (|A.34|) was guessed by comparing 
the expression ( ]63| ) for the screening mass in screened perturbation theory with the 
weak-coupling expression for m<, which is given in analytic form in Ref. 0. It was then 
verified numerically. 
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